Numerical evolution of radiative Robinson- Trautman spacetimes 



00 
O 

o 

(N 



H. P. de OliveirfEl and E. L. RodriguetQ 

Universidade do Estado do Rio de Janeiro 
Instituto de Fisica - Departamento de Fisica Teorica 
Cep 20550-013. Rio de Janeiro, RJ, Brazil 
(Dated: March 31, 2008) 

The evolution ol spheroids of matter emitting gravitational waves and null radiation field is stud- 
ied in the realm of radiative Robinson- Trautman spacetimes. The null radiation field is expected 
in realistic gravitational collapse, and can be either an incoherent superposition of waves of elec- 
tromagnetic, neutrino or massless scalar fields. We have constructed the initial data identified as 
representing the exterior spacetime of uniform and non-uniform spheroids of matter. By imposing 
that the radiation field is a decreasing function of the retarded time, the Schwarzschild solution is 
the asymptotic configuration after an intermediate Vaidya phase. The main consequence of the joint 
emission of gravitational waves and the null radiation field is the enhancement of the amplitude of 
the emitted gravitational waves. Another important issue we have touched is the mass extraction 
of the bounded configuration through the emission of both types of radiation. 
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I. INTRODUCTION 



An outstanding problem of classical general rel- 
ativity is the gravitational wave emission from the 
collapse of bounded distributions of matter to form 
black holes [l|. The knowledge of the underlying 
physical processes in which the nonlinearities of the 
field equations play a determinant role will be very 
useful in connection to the current efforts toward 
the detection of gravitational waves. Realistic mod- 
els of collapse must include rotation as well as the 
emission of radiation due to processes involving 
the matter fields inside the bounded source, there- 
fore constituting a formidable problem to be tackled 
numerically [SJ. In particular it will be of enormous 
importance to know the amount of mass radiated 
away by gravitational waves, the generated wave 
fronts and the influence of other radiative fields, such 
as electromagnetic and neutrino fields, on the emis- 
sion of gravitational waves. 

Recently, we have studied the full dynamics of 
Robinson- Trautman (RT) spacetimesjl] that can be 
interpreted as describing the exterior spacetime of 
a collapsing bounded source emitting gravitational 
waves The process stops when the stationary fi- 
nal configuration, the Schwarzschild black hole, is 
formed, and this happens for all smooth initial data. 
An interesting possibility in extending this model 
is to include a null radiation field that can repre- 
sent an incoherent superposition of electromagnetic 
waves, neutrino or massless scalar fields. There- 
fore, our objective here is to evolve the class of 
Robinson- Trautman geometries filled with pure radi- 



ation field^j that correspond to the exterior space- 
times of spheroids of matter which are astrophysi- 
cally reasonable distributions of matter. 

The paper is organized as follows. In Section 2 
we present the equations governing the evolution 
of radiative Robinson- Trautman spacetimes, along 
with the invariant characterization of gravitational 
waves through the Peeling theorem. In Section 3 we 
derive the initial data for Robinson- Trautman that 
represent the exterior spacetime of uniform and non- 
uniform spheroids of matter. The procedure of in- 
tegrating the set of nonlinear equations through the 
Galerkin method is outlined in Section 4. Although 
not presenting the details here, we have improved 
the implementation of the Galerkin method in which 
the issues of accuracy and convergence were dis- 
cussed in Ref. Q ■ Section 5 is devoted to the presen- 
tation and discussion of the dynamics of the model 
and the gravitational wave patterns. The issue of 
mass loss due to emission of gravitational waves is 
discussed in Section 6, where we suggest the valid- 
ity of a non-extensive relation between the fraction 
of mass extracted and the mass of the formed black 
hole. Finally, in Section 6 we present a summary 
of the main results as well as future perspectives. 
The linear approximation of the field equations is 
included in the appendix. Throughout the paper we 
use units such that SttG — c = 1. 



II. BASIC ASPECTS OF THE RADIATIVE 
ROBINSON-TRAUTMAN SPACETIMES 
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The Robinson- Trautman metric can be expressed 
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ds^ = (x(u,e)-^ + 2r—]du^ + 2dudr 
\ r K I 

-r^K^{u, e){de^ + sin^ ed(t)'^), (1) 

where m is a null coordinate such that u = constant 
denotes null hypersurfaces generated by the rays of 
the gravitational field, and that foliates the space- 
time globally; r is an affine parameter defined along 
the null geodesies determined by the vector d/dr 
and mo is a constant. The angular coordinates {9, (p) 
span the spacelike surfaces u — constant, r — con- 
stant commonly known as the 2-surfaces. In the 
above expression dot means derivative with respect 
to u and A(m, 9) is the Gaussian curvature of the 2- 
surfaces. Note that we are assuming axially symmet- 
ric spacetimes admitting the Killing vector d/d(t>. 

The Einstein equations for the Robinson- 
Trautman spacetimes filled with a pure radiation 
field are reduced to 



^("'^) = ]f2 --^ + ]^-]^cot0 (2) 



full nonlinear equation. The most general analy- 
sis on the existence and asymptotic behavior of the 
vacuum RT equation was given by Chrusciel and 
Singletonjl]. They have proved that for sufficiently 
smooth initial data the spacetime exists globally 
for positive retarded times, and converges asymp- 
totically to the Schwarzschild metric. Bicak and 
PerjesQ extended the Chrusciel-Singleton analysis 
considering the pure homogeneous radiation field, 
and established that for sufficiently smooth initial 
data the RT spacetime approaches asymptotically 
to the Vaidya solution 10] . In fact the spherically 
symmetric Vaidya metric is a particular solution of 
the field equations and Then, by assuming 
that K{u,9) = K{u), it follows from Eq. ^ that 
A(m, 9) = K~''^{u), and the RT equation yields 



K{u, 9) ~ K{u) cx exp 



6mo 



E'^{u)du 



The usual form of the Vaidya metric is recovered 
after a suitable change of coordinates, u ^ u = 
J du/K(u), for which the following expression for 
the mass function emerges 



-6mi 



k (Aesing)e 



(3) 



where the subscript 9 denotes derivative with respect 
to the angle 9. The function E'^{u, 9) appears in the 
energy-momentum tensor of the pure radiation field 
as 



T, 



fJ,l^ 



2 ^fl^L^f 



(4) 



where /p, = (5° is a null vector aligned along the de- 
generate principal null direction. The vacuum RT 
spacetime is recovered if E{u, 9) = 0, and the case 
in which the energy density does not depend on 9 
corresponds to the pure homogeneous radiation field. 

The structure of the field equations is typical of 
the characteristic evolution scheme: Eq. ([2]) is a 
constraint or a hypersurface equation that defines 
A(m, 9), and Eq. ([3]) or the RT equation governs the 
dynamics of the gravitational field. In other words, 
from the initial data K(uq,9) the constraint equa- 
tion determines X(uo,9), and after fixing the func- 
tion E'^{u,9), the RT allows to evolve the initial 
data. Therefore, this equation will be the basis of 
our analysis of the gravitational wave emission to- 
gether with a radiation field. 

Several important works have been devoted to 
the evolution of the vacuum RT spacetimes focus- 
ing on the issue of the existence of solutions for the 



Mvaidya{u) = 'mQK^{u). 



(5) 



By assuming further that the bounded collapsing 
matter emits radiation for some time until the 
Schwarzschild final configuration is settled down, its 
mass will be given by Mschw — "m-aKj, where Kf is 
the final value of K{u, 9). 

It will be important to exhibit a suitable expres- 
sion for the total mass-energy content of the RT 
spacetimes which we identify as the Bondi mass. To 
accomplish such task it is necessary to perform a 
coordinate transformation to a coordinate system in 
which the metric coefficients satisfy the Bondi-Sachs 
boundary conditions (notice that in the RT coordi- 
nate system the presence of the term 2rK /K does 
not fulfill the appropriate boundary conditions) . We 
basically generalize the procedure outlined by Fos- 
ter and Newman [lljl to treat the linearized problem, 
whose details can be found in Ref. O [ll. The 
result of interest is that Bondi's mass function can 
be written for any u as M{u,9) — moK^{u,9) + 
corrections, where the correction terms are propor- 
tional to the first and second Bondi-time derivatives 
of the news function; furthermore at the initial null 
surface u = uq these correction terms can be set to 
zero by properly eliminating an arbitrary function of 
9 appearing in the coordinate transformations from 
RT coordinates to Bondi's coordinates. Thus, the 
total Bondi mass of the system at u = uq is given 
by 
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M{uo) = -mo K^{uo,0) sm9d9. (6) 
2 Jo 

As a matter of fact, this expression will be very im- 
portant to our analysis of the mass extraction by 
gravitational waves. 

The invariant characterization of gravitational 
waves in RT spacetimes relies on the Peeling the- 
orem, which is based on the algebraic structure of 
the vacuum curvature tensor of the geometry ([!]). 
The Peeling theorem|U] states that the vacuum Rie- 
mann or the Weyl tensor of a gravitationally radia- 
tive spacetime generated by a bounded source ex- 
panded in powers of l/r has the general form 

n Nabcd , IIIabcd , Dabcd , 
Habcd ~ \ 1 ^3 H... {() 

where r is the parameter distance along the null vec- 
tor field d/dr, and the curvature tensor is evaluated 
in an appropriate tetrad basis f5|. The quantities 
Dabcd, IHabcd and Nabcd are of algebraic t ype 
D, type III and type N in the Petrov classification |16| 
and have the vector field d/dr as a principal null 
direction. The above equation contains the infor- 
mation about the radiation zone, that is, for large 
values of the distance parameter r the vacuum Rie- 
mann tensor is dominated by Nabcd /f, which is 
Petrov type TV. In other words, the field looks like 
a gravitational wave at large distances, in which the 
wave fronts are the it=const. surfaces with d/dr its 
propagation vector. In this vein, the non-vanishing 
of the scalars Nabcd is taken as an invariant cri- 
terion for the presence of gravitational waves in the 
wave zone. 

For the RT spacetimes filled with a pure radia- 
tion field the non- vanishing components of the Weyl 
tensor are 



W- 



2323 



W^0323 — 



-M^OlOl = 2Wo 



2too 



W^0303 — ^ 



D{u,0) A{u,e) + E^{u,e) 



(8) 



0202 



D{u,0) A{u,9)-E^{u, 



r 

where the functions 9) and D{u, 9) are given by 



Therefore the function D{u, 9) characterizes the 
wave zone and allows us to determine the evolution 
of the angular pattern of the emitted waves. It is 
important to remark that Gonna and Kramer [l3| 
have demonstrated, after performing a suitable co- 
ordinate transformation to the Bondi variables, that 
the news function of the RT metric is directly re- 
lated to D{u,9). Notice that both expressions are 
the same as found in the vacuum case, but the influ- 
ence of the radiation field on the gravitational wave 
amplitude is indirect since the evolution of K{u, 9) 
is altered by the inclusion of E{u, 9). 



III. INITIAL DATA: SPHEROIDS IN RT 
SPACETIMES 

A very important task is the construction of phys- 
ically plausible initial data K{uo, 9) for the RT equa- 
tion. We have chosen the family of prolate or oblate 
spheroids due to their astrophysical relevance, and 
also motivated by the works of Lin et al 17[ and 
Shapiro and Teukolsky[l^ that have analyzed the 
collapse of gas spheroids in Newtonian and relativis- 
tic situations, respectively. EardleyfT^ have focused 
on the collapse of marginally bound spheroids and 
the efficiency of the gravitational wave emission in 
this process. In order to construct the initial data 
correspondent to spheroids, wc follow the procedure 
outlined in Ref. [2^. We start with the oblate 
spheroid initial data whose initial step is to consider 
the oblate spheroidal coordinates {(, 9, (j)) defined by 



X ~ a-\/ + 1 sin 6* cos ( 
y = a\/(^'^ + 1 sin 9 sin (/ 



(11) 



a( cos ( 



where a is a constant, {x, y, z) are the usual cartesian 
coordinates and C^O, 0<6'<7r and < < 
2ti. Notice that the surfaces C ^constant represent 
oblate spheroids described by ^f^^^^^a) + ^jf^ — 1- 
The flat line element of the spatial surface E 



dx + dy + dz expressed in the 

Then, as 



above set of coordinates becomes ds^ — a 

cos2 9) ( + ^6*2) + (1 + C^) sin^ 9d(i 

the next step S is taken as the spacelike surface of 
initial data with geometry deflned by the following 
line element 



Aiu,9) 



1 



-As cot 1(9) 



K 



D{u,9) = 



1 d 
210 8^ 



^-^cot0-2f^ 
K K \ K 



(10) 



(3)(is2 = K^{C,9)\{C^ +COS 



(1 + C')sin2 



dC 



1+e 



d9' 



(12) 
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We assume time-symmetric initial data sucli that 
the vacuum Hamiltonian constraint [2 ij becomes 

V2'i' = 0, (13) 

where K — ^'^(C, 0) and is the flat space Laplace 
operator in oblate spheroidal coordinates. The sim- 
plest solution of this equation is 5'=constant that 
implies if=constant describing the Schwarzschild 
configuration. The general solution of this equation 
can expanded as a series of the product of Legendre 
polynomials of first and second kind[23|. By impos- 
ing suitable boundary conditions we obtain the fol- 
lowing non-fiat solution of that represents the 
exterior region of an oblate with uniform mass dis- 
tribution 

*(C,e) = + ^ («(C) + \mP2{cose))\ , (14) 

where a{() = arctan(l/C), /?(C) = (1 + 
3C^) arctan(l/^) — 3^, P2{cos9) is the second-order 
Legendre polynomial and Bq is a constant. The 
physical meaning of Bq becomes clear after rewrit- 
ing the metric (|12p with the above solution for large 
C (which implies in r a(), or 

^'^ds'^{l + ^+0{r--')ydsL,. (15) 

Therefore Bq is identified as the mass of the oblate 
spheroid as measured by an observer at infinity. 

At this point we can extract the initial data for 
the RT dynamics. Notice that the conformal func- 
tion evaluated at C = Co=constant contains all 
the information on the initial data correspondent to 
the exterior region of an oblate spheroid. Then, since 
the degrees of freedom of the vacuum gravitational 
field are contained in the conformal structure of 2- 
spheres embedded in a 3-spacelike surface, as estab- 
lished by D'Inverno and Stachel[23j that studied the 
initial formulation of characteristic surfaces, we are 
led to adopt the following initial data 



with foci at the points (0,0, ±a). After a straight- 
forward calculation following the same steps for the 
case of the oblate spheroids we obtain[2^ 

K{^o,e) = «[l + ^(«o - i^oP2(cose))]', (17) 

where ao = l/21n[(eo + 1)7(^0 - 1)] and /3o = 
l/2(3eo'-l)ln[(^o + 1)7(^0-1)] -3^0. 

It is possible to go further and establish a more 
general initial data that can be interpreted as de- 
scribing the exterior spacetime of a spheroid en- 
dowed with a non-uniform distribution of matter. 
It can be shown that if one considers a prolate 
spheroid, for instance, the initial data can be gener- 
ically expressed as 

anPn(cos6l)J , (18) 

where the closed expression for the additional term 
depends upon the way the matter inside the spheroid 
is distributed. When studying the consequences of 
such an initial data we will choose a suitable exact 
function of 9 instead the above generic expression. 

So far the initial data we have derived describe 
the vacuum exterior spacetime of bounded matter 
distributions. The derivation of the initial data 
for non-vacuum exterior spacetimes is a consider- 
able difficult task due to the fact that we have to 
solve nonlinear constraint equationsi21] instead of 
the Laplace equation whose analytical solution could 
be obtained. In this situation we assume that the 
same initial data as derived previously is extended 
to describe the exterior spacetime filled with pure 
radiation field, such that this model can account 
the collapse of a spheroid mass distribution emitting 
non-gravitational radiation. 

IV. NONLINEAR APPROACH: 
IMPLEMENTING THE GALERKIN 
METHOD 



^^(Co, e) = a[l + ^(^ao + i/3oP2(cos 61))] (16) 

with ao = a(Co) and /3o = PiCo)- 

A similar procedure can be done in order to derive 
the initial data representing the exterior spacetime 
of a uniform prolate spheroid. Now we start by in- 
troducing prolate spheroidal coordinates {£,,0, d)^^ 
X = a\/^2 — 1 sin cos 0, y = a-\/^^ — 1 sin sin (/>, 
z = aS^ cos 0, in which ^ > 1 and the surfaces 
^ — ^o=constant define prolate spheroidal spheroids 



In order to integrate numerically the field equa- 
tions ^ and ([3]) we have applied the Galerkin 
method [2^. For the sake of completeness we briefly 
outline here the procedure found in Refs. [5), [3] and 
improved in Ref. 0. Accordingly, the flrst step to 
establish the following Galerkin decomposition 

Kl{u,x) = A2e«("^-) = Agexp ( ^ hk{u)Pk{x)\ , 

\fe=0 / 

(19) 
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where a new variable x — cos 9 was introduced; N is 
the order of the truncation, and bk{u), k — 0,1, ...,N 
are the modal coefficients. The Legendre polynomi- 
als Pk (x) were chosen as the set of basis functions of 
the projective space, whose internal product is 



{Pj (x) , Pk (x)) = / P, {x)Pk (x) dx 



2 6. 



>kj 



2k + 1' 



The next step is to substitute Eq. lfT9|) into Eq. 
([2]) to obtain an approximate expression for X{u,x), 



\a{u, x) = 1 + E 2^(^ + l)bk{u)Pk{x) 

\ fc=0 / 

(20) 

Eqs. (fT9|) and (|20ll are now inserted into the RT 
equation ([3]) yielding the residual equation 



N 



Res(u,x) = Gtoq hk{u)Pk{x 

k=0 

[(l-x')X'J + 2E\u,x) 



1 

^0 



g-Q(«,a:) ^ 



(21) 



where the prime denotes derivative with respect to 
X. In general the residual equation approaches to 
zero when N 00 meaning that the decomposi- 
tion tends to the exact K{u,x). The Galerkin 
method establishes that the projection of the resid- 
ual equation onto each basis function Pn{x) vanishes, 
or {Res{u, x), Pn{x)) = 0. These relations can be 
cast in the following way 



-S^^(E'{u.^).P„{a:)), (22) 

where n = 0,1,2.., N. Therefore, the dynamics of 
RT spacetimes is reduced to a system of {N + 1) 
ordinary differential equations for the modal coeffi- 
cients, or simply a dynamical system. It is worth of 
mentioning that the issues of accuracy and efficiency 
of the dynamical system approach are discussed in 
Ref. 0. 

To integrate the dynamical system ((22|) we need to 
provide the initial conditions bk{uo), k = 0,1, 2, .., N 
that represent the initial data K[uq,x) constructed 
in the last Section (Eqs. ^ and ^). These 

initial values are obtained from the Galerkin decom- 
position of _fi'(uo, x) (cf. Eq. (19)), which yields 



2{liiK{uo,x),Pj] 



(23) 



As we have mentioned in the last Section the sim- 
plest form for the initial data correspond to the 
Schwarzschild solution K{uq, x) — constant that is 
described exactly using the above decomposition in 
which the only non- vanishing modal coefficient is 69 • 
In order to illustrate the effect of increasing the trun- 
cation order N we have exhibited the plots of the 
relative error, R.E., between the approximate and 
exact expressions for the prolate spheroid described 
by Eq. that is 



R.E. 



K{uo,x) — Ka{ui),x) 



K{uo,x) 



Note that even for a small truncation order TV = 4 
the error is acceptable. In the numerical experiments 
of the next Section we have fixed N — 14. 
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FIG. 1: Plot of the relative error between the approxi- 
mate and exact initial data corresponding to a prolate 
spheroid for a — l,^o ~ 1-01 and Bo = 2. The curves 
correspond to truncations orders A'' = 4, 8, 14, 20. Notice 
that even for the modest truncation A'^ = 4 the error is 
acceptable. 



V. NUMERICAL RESULTS 

As we have shown to integrate the set of equa- 
tion it is necessary to specify the initial con- 
ditions determined by the initial data K(uo,x) (cf. 
Eq. [^5)1 . and the function E(u,x). We start with 
uniform spheroid initial data emitting homogeneous 
radiation (non-gravitational). It is reasonable to 
impose that E{u) is an arbitrary decreasing func- 
tion of the retarded time u. Amongst several pos- 
sible choices [2^ for E{u) we may select E^{u) — 
Equ" exp(— /3u), where Eq, a and /3 are positive con- 
stants; in particular if a ^ we describe a pulsed 
radiation. For any function E{u) it is important to 
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FIG. 2: Evolution the modal coefficients bo{u) and &4(«) 
for the prolate spheroid with a — 1, Bo = 0.1 and 
^0 ~ 1.01 for the vacuum and non-vacuum cases. We 
have considered the homogeneous radiative field E'^{u) — 
Egue~°'^^ , where for the first and second figures (a) 
Ei = 0.05, 2.0 and (b) E^ = 0.1, 5.0, respectively. 

remark that the influence of the radiation field is ba- 
sically due to the term (2£'2(u), P„(a;)) = 4E'^{u)Sno 
(see Eq. (P^ ) that only appears explicitly in the 
equation for b(){u)] the influence on other modal co- 
efficients is indirect through the nonlinear couplings 
between bQ{u) and the remaining ones. 

Let us consider the initial conditions 6fc(0), k = 
0, 1, .., iV that correspond to the prolate spheroid de- 
scribed by Eq. (fTT]) . By inspecting this expression 
only even Legendre polynomials are present, and 
therefore only the even modal coefficients are ini- 
tially distinct form zero. According to our previous 
work on the dynamics of vacuum RT spacetimesjH, 
we have shown that the odd modal coefficients re- 
main unaltered while the even modal coefficients 
tend to zero asymptotically. The exception is Bq^u) 
that approaches to a constant value, consequently 
characterizing the Schwarzschild solution as the end 
state, which in agreement with the analytical works 
of Chrushiel and Singleton Q. The Schwarzschild 
mass of the final configuration is given by 

Mbh ^ moAlei'^^^f^ (24) 

where bo{uf) is the asymptotic value of the zeroth 
modal coefficient, and u/ is large but finitc[36^. The 
final mass Mbh is always smaller than the initial 
mass M{uo) as should be expected since part of the 



initial mass is carried away by gravitational waves. 
We will return to this issue in Section 6. 

We have performed several numerical experiments 
using the set of equation (I22|) with truncation or- 
der N = 14, and homogeneous radiation field E'^{u) 
with distinct values of E^o , ot and (3. Physically, we 
have introduced a non-gravitational way of extract- 
ing mass that was discussed for the first time by 
Stachel[27|. but it can be also seem from the lin- 
earized formula of mass loss (jA5l) . In this way we 
expect that the mass of the final configuration be 
smaller than the similar one formed in the absence 
of the homogeneous radiation field. As a matter 
of fact, the numerical results have confirmed this 
feature; the Schwarzschild solution is the asymp- 
totic configuration, that is b2k 0, k > 1, &o ^ 
bo{uf) = constant and smaller than found in the 
vacuum case. Another important consequence of the 
homogeneous radiation field is that the approach to 
the Schwarzschild solution occurs faster than in the 
vacuum case. In Fig. 2 we have illustrated these fea- 
tures with the plots of the modal coefficients 6o(u) 
and 64 (u). 




^-2 024 -a. n o d 



FIG. 3: Dynamics of all modal coefficients through the 
two-dimensional plots of K{u, x) starting from the pro- 
late spheroid initial data (a = 1, Bo = 0.5, .^o ~ 1.05) at 
uo — 0, and for several times until the Schwarzschild fi- 
nal configuration (represented by a circle) is formed. The 
first graph refers to the vacuum RT spacetimes, while 
the second includes the homogeneous radiation field 
E'^{u) = 4.0exp(— 0.5u). One of the effects of the ho- 
mogeneous radiation field is to produce a Schwarzschild 
configuration with smaller mass (smaller circle). 

The dynamics of all modes together is shown 
in Fig. 3 for the vacuum and non- vacuum cases. 
Both graphs consist in sequences of polar plots 

of Kaiu,x) = Aoexp (l/2J2k=o ^kiu) Pkix)j ^^"^ 
times u until the Schwarzschild configuration (rep- 
resented by a circle) is settled down. As expected, 
the presence of radiation produces a smaller circle 
due to the additional mass extracted. 

The emission of gravitational waves consti- 
tutes the central physical aspect of the radiative 



7 



Robinson- Trautman spacetimes. As we have seen 
the invariant characterization of the presence of 
gravitational waves was estabhshed on the basis of 
the Peehng theorem, where D{u,9) given by Eq. 

contains aU information about the angular and 
time dependence of the gravitational wave amplitude 
in the wave zone. This function can be conveniently 
expressed in terms of a; = cos 6 and Q{u, x) by 




FIG. 4: Polar plots showing the evolution of the angular 
pattern of gravitational waves emitted at the radiation 
zone for the vacuum (left panel) and non-vacuum cases 
(right panel), where the quadrupole mode is dominant. 
The arrows indicate the first instants denoted by no = 
0, ui = 0.06, U2 = 0.18 in which the angular pattern 
of D{u,x) were evaluated. The homogeneous radiation 
field described by E^{u) = exp(-0.5ii) with = 20 
does not afi^ect the pattern of gravitational waves but 
increases the peaks as indicated above. In both cases 
the initial conditions correspond to a prolate spheroid 
with a = 1, ,^0 = 1-01, Bo = 0.1. 

The angular distribution -D(m, x) at each instant 
u is determined from the Galerkin decomposition 
Q{u,x) — X^fcLi bk{u)Pk{x) along with the dynam- 
ical system for the modal coefficients. Although not 
appearing explicitly in the above expression, the ra- 
diation field alters the behavior of the modal coeffi- 
cients and consequently the function D(u,x). 

A very useful way of depicting pattern of the grav- 
itational waves at the radiation zone is through by 
polar plots of D{u,x). In Fig. 4 we present such 
polar plots evaluated at the initial stages of evolu- 
tion (as indicated by the arrows) and correspond- 
ing to the prolate spheroid initial data p?)) for the 
cases of vacuum and homogeneous radiation field 
E^{u) = i?o exp(— 0.5m). The patterns of emit- 
ted waves have symmetry of reflection by the plane 
6 = tt/2 {x = 0), which is a direct consequence of 



= tt/2 
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FIG. 5: Polar plots showing the evolution of the angu- 
lar pattern of gravitational waves emitted by an oblate 
spheroid (|16p for vacuum (upper panel) and non-vacuum 
cases (lower panel). The structure now resembles of 
that octupole pattern, and as in the case of the prolate 
spheroid the homogeneous radiation field increases the 
peaks of gravitation waves. Again the arrows indicates 
distinct instants Mo = 0,ui = 0.3, U2 = 0.6. 

the absence of all odd modes 62fe+i- Indeed, the 
symmetry of the pattern is consistent with the final 
configuration being a Schwarzschild solution at rest 
with respect to a distant observer, since no net mo- 
mentum is carried out by the gravitational waves. In 
Fig. 5 similar plots are depicted taken into account 
oblate spheroids initial data. 

Two interesting aspects are worth of remarking. 
The first is the predominant quadrupole structure of 
the pattern of gravitational waves, and the second is 
the influence of the homogeneous radiation fleld that 
does not alter the pattern of gravitational waves but 
enhance the amplitude of the peaks. To explain this 
result we might remember that the radiation field 
induces a faster decay of the modal coefficients (cf. 
Fig. 2), and consequently produces higher values of 
bk- From Eq. (^5]) . D{u,x) depends directly on the 
derivatives of the modal coefficients, and therefore 
its value will be eventually increased. In order to 
show quantitatively the amplification of the gravi- 
tational wave amplitude due to the radiation field, 
we have evaluated the total flux of gravitational ra- 
diation in the wave zone at each instant, which is 
proportional to 

a{u) = J \D{u,x)\dx, (26) 

according to a distant observer. In Fig. 6 (y{u) is 
plotted for the vacuum and non-vacuum RT space- 
times indicated by the labels (a), (b) and (c) for 
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increasing values of Eq taking into account that 




250 300 



FIG. 6: Influence of the radiation field E^{u) = 
-Bq'^ exp(— 0.5u) on the behavior of the flux of emitted 
gravitational waves. We have selected EI = 3.0, 4.0, 5.0 
in (a), (b) and (c), respectively. The initial data is the 
prolate spheroid ([T7)) with a = l,^o = 1.01, Bo = 1.0. 



We now turn our attention to a more general ini- 
tial data that represent the exterior spacetime of a 
non-uniform prolate spheroid given by Eq. (jlSp . 
Before discussing the numerical results, it is intu- 
itive to expect that the collapse will proceed in a 
non-uniform way, meaning that distinct parts of the 
spheroid will collapse faster/slower than others. As 
a consequence, the pattern of gravitational waves 
should not display the same symmetry of reflection 
by the plane 6 — tt/2 {x — Q) produced in the case 
of an uniform spheroid (cf. Figs. 4 and 5). In this 
case there will be a net flux of momentum carried 
out by gravitational waves that would be responsible 
for the recoil of the spheroid. 

In order to study more quantitatively this feature, 
we consider, for the sake of convenience, that the last 
expression present on the rhs of Eq. (|18p is writ- 
ten as X^^^o <^nPn{x) = Aoexp(-(a; - xq)^), and 
that reflects somehow the description of the non- 
uniformity of the distribution of matter inside the 
spheroid, or alternatively, a perturbation of the uni- 
form distribution. As a consequence, the Galerkin 
decomposition yields that all modal coefficients bk (0) 
are initially distinct from zero. We have first evolved 
the dynamical system for the vacuum RT space- 
times {E(u,x) = 0), and contrary to what we have 
found previously, all modal coefficients tend to con- 
stant values asymptotically (for large u) as illus- 
trated in Fig. 7(a) with the plot of the modal co- 
efficient hx{u). In Fig. 7(b) we have depicted some 
of the numerically resulting configurations Ka{x) = 

exp^l/2X^Q bk{u f) Pk{x)^ . In order to identify 

these stationary configurations we recall that the 
other possible static solution of the RT equation is 



-1 ^[15 ^ 05 1 

X 

(b) 

3.4e-08 - 
3.2e-08 - 
3e-08- + 
t 2.8e-08 - 
2.6e-08 - 
2.4e-08 - 
2.2e-08 - 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
(c) 

FIG. 7: (a) Behavior of b\{u) indicating that its asymp- 
totic value is distinct from zero. The initial data are non- 
uniform oblate spheroids characterized by i^o = l,Bo = 
1 and Y^'^^Q a„P„{x) = Aoexp{—{x — 0.3)^), where 
{i)Ao = 0.7, lii)Ao = 0.2, {iii)Ao = 0.05. (b) Plots of the 
final configuration Ka{x) = exp{l/2'}2 bk{uf)Pk{x)). 
These curves are quite well fitted by the analytical ex- 
pression (|27[) after choosing appropriate values for Kq 
and fi. In (c) we present the rms error between Ka and 
K{x) for each value of Aq. 



cosh n + X sinh fj. ' 

where Kq and /i are constants. This solution was 
found by Bondi et al[lH| and represents a boosted 
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black hole traveling with velocity v ~ tanh/i along 
the symmetry axes. The total mass of this configu- 
ration is 

Mq = mrest cosh^ = J"*'""* , (28) 
V 1 — w 

where mrest = ttioKq is the rest mass = 0) with 
respect to a distant observer. Fig. 7(c) shows the 
rms error between each numerical solution of Fig. 
7(b) and the exact solution P7|) after convenient 
choices of Kq and fi. Therefore, if all modes are ini- 
tially excited the asymptotic solution is a boosted 
black hole. 

It is interesting to notice that in the present sit- 
uation the bounded source is accelerated from the 
rest at u = until the emission of gravitational 
waves finishes. According to Fig. 8, the pattern 
of gravitational waves exhibits two prominent lobes 
whose symmetry line is the direction of recoil; also 
the inclination of the lobes increases with time (cf. 
Fig. 8)) together with the decrease of their am- 
plitudes. These features are typical of the pattern 
of waves emitted by an accelerated point particle 
and are known to form a pattern of bremsstrahlung 
radiation [Sm. 




FIG. 8: The graph on the left shows the polar plot of 
D{u, x) exhibiting a typical bremsstrahlung pattern of 
the gravitational radiation, whereas the detail of one of 
the lobes evaluated nt u — 0, u = 0.06 and u = 0.18 is 
depicted on the right. Note that two effects are present: 
the amplitude of the gravitational waves decreases, and 
the inclination of the lobe changes forming a slightly less 
sharp forward cone in the direction of motion, as ex- 
pected from the usual bremsstrahlung radiation. 

The final task is to consider inhomogeneous ra- 
diation fields which seems to be more consistent 
in the realm of non-uniform spheroids. Although 
an inhomogeneous radiation field might arise as a 
consequence of general internal processes inside the 



spheroid where the non-uniform matter distribution 
may play a relevant role, we have not found any 
physically realistic function to describe such form 
of radiation. Nonetheless, we can make some con- 
siderations about this situation. The dynamics of 
the modal coefficients will be directly affected by 
an inhomogeneous radiation field, as it can be seen 
from the projection ^i5^(it, x), P„(x)) that in gen- 
eral does not vanish for any n. Possibly, the decay of 
the modal coefficients to their final values will done 
faster than in the absence of radiation, as in the 
case of homogeneous radiation. Besides the addi- 
tional mechanism of extracting mass of the spheroid 
provided by the radiation field, there will be also a 
additional contribution to the total balance of lin- 
ear momentum. Therefore, the final configuration 
will be in general a boosted black hole with the fol- 
lowing characteristics: its mass is smaller and its 
velocity can be greater or smaller, respectively, than 
the mass and velocity of the corresponding boosted 
black hole found without radiation. Concerning the 
pattern of gravitational waves, we expect that con- 
trary to the homogeneous case, the inhomogeneous 
radiation field will change the produced pattern of 
waves. 



VI. THE MASS LOSS 

We have already studied the issue of mass loss 
due to the emission of gravitational waves before 0, 
[T^ . [20| in the realm of vacuum Robinson- Trautman 
spacetimes. The problem consisted in determining 
numerically the relation between the fraction of mass 
lost A, or the efficiency of the process, with the final 
mass of the configuration Mbh- Here A is defined 

_ \M{uo)-Mbh\ 
^ = MM ' ^ ^ 

where Af (mq) is the initial mass evaluated according 
to Eq. ([6]). As the main result we have shown that 
the set of numerically generated points (Mbh,^) 
can described by the following expression inspired 
from the Tsallis statistics[31], 

A = A„ax(l-2/"^)T^, (30) 

where y — Mbh / fno^Q, Amax is the maximum ef- 
ficiency, 7 and q are the other parameters of the 
distribution law. We have found that for distinct 
initial data the best fit required that g sa 0.5, the val- 
ues of 7 and Amax depend on the particular choice 
of the initial data. Here the idea is to revise these 
results using our improved numerical method 7], to 
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present some analytical work concerning the relation 
between Mbh and A, and to study the consequences 
of introducing a radiation field. 

Before presenting the analytical and numerical re- 
sults it is important to stress that Eq. ([50]) may not 
be considered an artificial relation connecting A and 
Mbh- To support this point of view we recall that 
a black hole can be understood as thermodynami- 
cal system in equilibrium, therefore it is a state of 
maximum entropy. The question we address is the 
following: given an initial data K{uo, x) now viewed 
as a "non-equilibrium thermodynamical state", how 
much mass must be extracted through gravitational 
waves such that "thermodynamical equilibrium" or 
the black hole formation is achieved? The empirical 
answer to this question is given by Eq. ([50)) . Never- 
theless, the choice of this function is inspired in two 
aspects: non-extensive relations involving thermo- 
dynamical quantities are typical of systems in which 
a long-range interaction such as gravitation is the 
relevant onefs^], and the black hole entropy is a non- 
extensive quantity [1^. 

Let us consider the regime of small mass extrac- 
tion or equivalently A <^ Amax- As a matter of fact 
this regime can be approximately described by the 
linearized equations (cf. Appendix A). Therefore, by 
setting E'^{u,x) = into the equation that governs 
the loss of mass (Eq. (jASp ). A can be calculated 
analytically: 

k 

where e(w/) is the final value of eo(w)- Notice 
that A ~ C(e^) implying that Seo = eo{uf) — 
eo(0) ~ C(e^) (indeed, the linear approximation 
dictates that Sea = 0). On the other hand, y — 
MBH/maAl « 1 + 3/2eo(w/) - 1 + 0(e). Combin- 
ing these results it follows that 

y- 1 cx Ai/2 ^ y w 1 + const.Ai/2. (32) 

By noticing that the nonextensive relation (|30p ren- 
ders y « l-|-(A/A,„ax)^~V7 in the limit A < A^ax, 
we conclude that q = 0.5, which is reproduced by 
the numerical result. Therefore, no matter is the 
initial data the parameter the regime of very small 
efficiency dictates that q = 0.5. 

In Fig. IHa) we have depicted the plots of the 
points {Mbh,^) generated numerically with the 
corresponding best fits of the function ([50]) (con- 
tinuous lines). The oblate spheroid was taken 
as the initial data where -Bo is the free parame- 
ter while a = 1 and C = 0.01,0.5 and 1.0. Note 
that C is associated to the oblateness through / = 
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FIG. 9: (a) Log-linear plots of the best fit between the 
points obtained numerically and the non-extensive func- 
tion (|3Up . The initial data correspond to uniform oblate 
spheroids with a — 1, = 1.0 (circle), ( = 0.5 (di- 
amond), = 0.01 (box), and Bo is the free parame- 
ter. The best fit required q ~ 0.5003, q ~ 0.5007 and 
q ~ 0.5024, respectively, (b) Plots of the fraction of 
mass loss considering uniform prolate spheroids as the 
initial data whose exterior is vacuum RT and homo- 
geneous radiation field E^{u) = Eq exp{—0.5u) for (a) 
El = 0.001, (h) El = 0.01 and (c) El = 0.05. In both 
cases y — Mbh /moAo. 



1 ~ C/\/l + C^i a measure of how much the ini- 
tial distribution of matter deviates from the sphere. 
Therefore, according to Fig. Ufa) the maximum ef- 
ficiency attained increases with the oblateness (or 
decreases with () as should be expected. The nu- 
merical values obtained for the maximum efficiency 
are A^ax ^ 0.01038, 0.0510, 0.2417, respectively for 
C = 1.0,0.5,0.01. It is worth noting that the last 
maximum efficiency is quite high due to the abnor- 
mally enlarged degree of oblateness. For sake of com- 
paration, in the realm of our Solar System the high- 
est oblateness belong to Saturn with /satmn — 0.097 
(C ^ 2.08), and Jupiter /j^pitcr ^ 0.065 (C ^ 2.64), 
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which are quite smaU. 

We have noticed that one of the most important 
consequences of the radiation field is to enhance 
the emission of gravitational waves, increasing, as 
a consequence the mass loss. The numerical rela- 
tion between Mbh and A is shown in Fig. [HJb), 
where we have considered for the homogeneous ra- 
diation field £"^{11) = E^exp{-0.5u). The effect of 
increasing Eq is indicated by the successive devi- 
ations of the points corresponding to the vacuum 
case. Although we could not find a function to fit 
the data, it is possible to estimate analytically the 
fraction of mass extracted in the late stages domi- 
nated by the radiation field, which is described by 
the Vaidya solution. For the sake of simplicity we 
consider that asymptotically efe — > 0, A; > 1 and 
eo — *■ constant. The late time behavior of eo{u) 
can be described approximately by the linear ap- 
proach such that Co ~ — 1/3too-E^(u), which yields, 
Scq w — i?o/3mo/3. It can be shown that the fraction 
of mass loss is A w 1 — exp{—EQ/2moP), which is 
in agreement with the numerical data as shown in 
Table 1. 
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Relative Error 


1 


0.0951625820 
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0.0009995002 
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0.005 
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0.001 


0.0000999950 
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VII. FINAL REMARKS 

In this paper we have studied the dynamics of 
the exterior spacetime of a bounded source emit- 
ting gravitational waves and a null radiation field 
characterized by the function E{u,6) (cf. Eq. 2]). 
The joint emission gravitational and some other type 
of radiation must be expected in a realistic non- 
spherical collapse. In order to provide a simple 
model with such features, we have derived the ini- 
tial data that represent the exterior spacetimes of 
uniform and non-uniform spheroids of matter in the 
realm of Robinson- Trautman geometries. 

Basically, the main implication of introducing the 
null radiation field is to produce a faster approach to 
the final configuration if compared with the case of 
vacuum {E{u, 6) =0), provided that E{u, 9) is regu- 
lar and decays with the retarded time u. Indeed, as 
a consequence of this rapid relaxation towards the 
stationary configuration, the amplitude of the emit- 
ted gravitational waves is enhanced as illustrated in 
Figs. 4 and 5. In Fig. 6 a more quantitative view 



of this effect is exhibited. A similar feature can be 
found, for instance, in the process of formation of su- 
pernova in which asymmetric distribution of neutri- 
nos can pro duce a growth in the gravitational wave 
emission [3Jj. 

The angular distribution of gravitational radia- 
tion at the wave zone bestows important information 
about the initial data and the final configuration. In 
this vein, distinct initial data will produce distinct 
pattern of emitted waves, and this fact is clear from 
Fig. 4 and 5 whose initial data arc uniform prolate 
and oblate spheroids; also the pattern of Fig. 8 is 
due to non-uniform prolate spheroid. On the other 
hand, if the pattern has symmetry of reflection by 
the plane = 7r/2 the final configuration will be 
the Schwarzschild black hole, but if such symme- 
try is not present as shown in Fig. 8, a boosted 
black hole will be the end state. As we have seen, 
the emission of gravitational waves transport a net 
amount of linear momentum, which is responsible 
for the recoil of the isolated source due to energy- 
momentum conservation. In this process the mat- 
ter source is accelerated producing typical pattern of 
bremsstrahlung |30| radiation of gravitational waves. 
The recoil of the remnant (a black hole or a pulsar) 
resulting from a non-spherical collapse is expected, 
and has been studied largely in the literature, mainly 
in the context of supernova explosions [ssj. In par- 
ticular, according to Eardley[23| pure gravitational 
wave recoil in the collapse applies only to the case 
of collapse to a black hole. 

We have also revised the important issue of mass 
extraction due to the emission of gravitational waves 
in the realm of vacuum RT spacetimes. Our im- 
proved code has confirmed our previous numerical 
result concerning the non-extensive relation between 
the fraction of mass extracted and the black hole 
mass. As a new result we have derived analytically 
that q = 0.5, which was confirmed by the numerical 
experiments. 

Finally, the next step of our investigation will 
be to consider the dynamics of more general ax- 
isymmetric spacetimes described in the Bondi prob- 
lem, where a numerical code based on the Galerkin 
method is under development. One of our interests 
is the mass loss due to gravitational wave extraction. 

The authors acknowledge the financial support of 
the Brazilian agencies CNPq and CAPES. We are 
also grateful to Dr. I. D. Soares for useful comments. 

APPENDIX A: EXACT SOLUTION OF THE 
FIELD EQUATIONS: LINEAR 
APPROXIMATION 

It is useful to study the linear regime of the field 
equations of RT spacetimes endowed with radiation 
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field, which is a direct generalization of the Foster 
and Newman[ll[ solution for the vacuum case. Basi- 
cally, we follow the approach of Foster and Newman 
by assuming that K(u, 9) is approximated as 



A' 



n=0 



i{u) P„(cos 0) 



, (Al) 



where Aq is a constant, |e„| ^ 1 and P„(cos 9) is 
the Legendre polynomial of order n. For the vacuum 
case Foster and Newman have chosen n > 2 with- 
out loss of generality, since n — Q and n = \ both 
represent spheres (considering, of course, the sub- 
manifold K"^ {d9'^ + sin^ 9 dc/)'^)) that in the latter 
case has its center displaced at a distance ei Aq along 
the axis of revolution. The linearized field equations 
are reduced to 



A:=0 



En Pn{cOS 6*) = - ^ fc„e„ P„(cOS 



3mo 
(A2) 

where fc„ — -jj^ (n -|- 2)!/(n — 2)! is positive for all 
n > 2. Let us consider first the homogeneous radi- 
ation field, E^{u,9) = E'^{u), for which the above 
equation can be rearranged as 



eb(w) 



3mo 



£'^(m), e„{u) = -fc„e„(u). (A3) 



The second equation renders immediately e„(u) = 
e„(0) exp(— /c„m), e„(0) being the initial value of e„ 
evaluated on the hypersurface u — uq — 0. From 
this analytical solution all coefficients e„ — > as 
u — > oo, which is the same outcome of the vac- 
uum RT geometries. The zeroth mode behaves as 
eo{u) = eo(0) - /"^ E'^{u)du, indicating the ap- 
proach to the Vaidya spacetime asymptotically. In 
this case, the corresponding line element that char- 
acterizes the Vaidya geometry can be obtained af- 
ter the same coordinate transformation indicated in 
Section 2. 

By considering an inhomogeneous radiation field, 
we may set E^{u,9) — J2jLo fj{'^)Pj{'^'^^(^) to 
model a physically reasonable form for the energy 
density of the radiation. Taking into account this 
expression into Eq. (jA2[) . the equations for each 
mode e„(w) can be integrated yielding 



,{u) = e„(0)e" 



-knU 



-knU 



3mn 



/„(u)e'=""dM. 

(A4) 

If n = the solution corresponding to the homo- 
geneous radiation case provided E'^[u) — fo{u) is 
recovered. Again, we also may set fi(u) — with- 
out loss of generality. For those modes n > 2, the 
first term on the rhs vanishes in the limit u ^ oo, 
and the long term behavior of e„(w) is governed 
by the second term. It is possible to know the 
outcome of this asymptotic limit without specify- 
ing the functions /n(u), n > 2. For instance, sup- 
pose that all fn{u) are regular in some interval of 
time say, uq < u < ui, representing a pulsed emis- 
sion of radiation, or even be decreasing functions 
such that lim.u^oo fn{u) — 0. As a consequence, 
the second term on the rhs vanishes asymptotically 
rendering again e ^ as m ^ cxd. Under these 
conditions and, in addition assuming further that 
lim„_,oo E^(u)du is finite, the asymptotic state 
will be a Schwarzschild black hole [37]. 

As a final piece of our study we generalize the orig- 
inal result of Foster and Newman on the change in 
the total mass-energy of the system. We have fol- 
lowed an alternative approach (instead of using the 
Newman-Penrose formalism), which consists in ex- 
pressing the RT equation in terms of the mass func- 
tion M{u) ~ l/2mo /J" K^{u,9)sm9d9, and taking 
into account the linearized solution (Eqs. (|Aip and 
(|A2[) '). After some steps and considering the homo- 
geneous radiation field for the sake of simplicity, we 
have obtained 



dM{u) 
du 



IQAq 



E 



(" + 2)! el{u) 



in 



A3 



2)! 2n 
3 



1 



^E^{u)il + -eo{u) 



(A5) 



The first term on the rhs, found by Foster and 
Newman, means that the mass loss by gravitational 
wave extraction is a second order effect, while the 
second term accounts for the mass change due to 
the radiation field. In this sense, we can interpret 
the present linearized solution as representing as a 
bounded source emitting gravitational waves and a 
null radiation field. 
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